Thermally activated Peierls dimerization in ferromagnetic spin chains 
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We demonstrate that a Peierls dimerization can occur in ferromagnetic spin chains activated by 
thermal fluctuations. The dimer order parameter and entanglement measures are studied as func- 
tions of the modulation of the magnetic exchange interaction and temperature, using a spin-wave 
theory and the density-matrix renormalization group. We discuss the case where a periodic modu- 
lation is caused by spin-phonon coupling and the case where electronic states effectively induce such 
a modulation. The importance of the latter for a number of transition metal oxides is highlighted. 
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Structural instabilities of electronic systems can occur 
due to the coupling of electronic and lattice degrees of 
freedom (phonons). They are particularly important for 
quasi one-dimensional (ID) systems where the gain in 
electronic energy due to a lattice distortion often out- 
weighs the cost in elastic energy. A well known example 
is the Peierls instability [l| of the ID free electron sys- 
tem towards a static lattice distortion determined by the 
Fermi momentum. For a commensurate distortion, an 
excitation gap is opened turning a metallic system into 
a band insulator. This Peierls metal-insulator transition 
plays an important role, for example, in organic charge- 
transfer solids [2]. A related instability occurs for an- 
tiferromagnetic (AFM) spin chains coupled to phonons. 
Here magnetic energy is gained by distorting the lattice 
which can lead to the so called spin-Peierls (SP) transi- 
tion. Although a SP phase transition was first observed 
in organic materials [3(, it was the discovery of such a 
transition in CuGeC>3 by Hase et al. [4[ that has led to 
great interest in these phenomena [5| . 

Quite recently, another type of Peierls instability for 
spin chains has been found which is not driven by spin- 
phonon coupling but rather by a coupling of the spins 
with electronic degrees of freedom (orbitals) . Here a fer- 
romagnetic (FM) spin chain shows a periodic modulation 
(dimerization) of the magnetic exchange in a certain fi- 
nite temperature region while the ground state is the 
uniform fully polarized FM state [6( . In [7], |8| it has been 
argued that this mechanism is responsible for the remark- 
able physical properties of YVO3 in the finite tempera- 
ture C-type AFM phase. Clearly, the lattice will react to 
a modulation of the magnetic exchange, however, spin- 
phonon coupling is not the driving force in this case and 
any lattice distortion is only a secondary effect. 

In this Letter we want to establish general mechanisms 
which can drive a Peierls dimerization in FM spin chains. 
To highlight the differences between AFM and FM chains 
we will first consider a coupling to lattice degrees of free- 
dom. The phonons are often treated adiabatically which 
is justified if the phonon frequency is smaller than the 
Peierls gap. In the adiabatic approximation the Hamil- 



tonian can be written as H 
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and E c \ — NK5 2 /2. Here J is the exchange constant, Sj 
is a spin S operator at site j, and N is the number of sites. 
In the absence of a magnetic field the modulation is ex- 
pected to be commensurate with wave vector k = n and 
is parameterized by S e [0, 1]. The elastic energy E c \ con- 
tains the elastic constant K. Both parameters K, 6 can 
be related to the spin-phonon interaction strength. Note 
that writing the Hamiltonian as a sum of a magnetic and 
an elastic part as in Eq. ([I]) corresponds to the random- 
phase approximation (RPA) by Cross and Fisher [9( . Al- 
though the model (JTJ is strictly ID, the static, mean-field 
(MF) treatment of the three-dimensional phonons allows 
for a finite temperature phase transition if S(T) is treated 
as a thermodynamical degree of freedom determined by 
minimizing the free energy. 

Let us start with the case where Sj-Sj+i — » SjSJ +1 + 
S]S? +1 = (S+S- +1 + SjS+ +l )/2, i.e., we replace the 
SU(2)-symmetric spin exchange by an XX-type of in- 
teraction. In this case the sign of J does not matter and 
the system becomes equivalent to a free spinless fermion 
model by Jordan- Wigner transformation. The Hamilto- 
nian is then easily diagonalized by Fourier transformation 
and in the ground state for small S one finds a gain in 
magnetic energy E mag ~ S 2 \nS. This outweighs the cost 
in elastic energy E c \ ~ 8 2 and constitutes the Peierls 
instability for lattice fermions [id] . For the isotropic 
antiferromagnet (J > 0, SU(2)-symmetric exchange), 



field theoretical arguments show that E n 



.J4/3 



Again this outweighs the cost in elastic energy leading to 
a SP transition and the opening of a spin gap, A ~ <5 2 / 3 . 
Contrary to the two cases discussed above there is no 
gain in magnetic energy in the ground state for FM cou- 
pling, J < 0. For 5 G [0, 1) the ground state is always the 
fully polarized FM state. We will show in the following 
that thermal fluctuations can, however, activate a Peierls 
dimerization. We will use the density-matrix renormal- 



2.5 



1.5 



0.5 




0.1 \ K 0.2 



K/IJI 



0.3 



0.4 



0.5 



FIG. 1: (color online) Phase diagrams for the dimerized AFM 
and FM S = 1/2 Heisenberg chains. The dot-dashed line 
depicts the second order SP transition from the uniform (U) 
to the dimerized (D) phase for the AFM chain. For the FM 
chain, the D phase exists only at finite T (in units of fcs = 1) 
and only if K < K c ~ 0.118 \J\ - — the transition is either 
second or first order, as shown by solid and dashed lines, 
respectively, and changes its character at the tricritical point 
(TCP). The insets show the order parameter S(T) for: (a) 
the FM chain with K/\J\ - 0.1, and (b) the AFM chain with 
K/J — 2 (the lines are guides to the eye). 



ization group applied to transfer matrices (TMRG) to 
study this effect. The TMRG algorithm is based on a 
mapping of the ID quantum onto a two-dimensional clas- 
sical system. A transfer matrix is then defined allowing it 
to perform the thermodynamic limit exactly, i.e., all the 
numerical results presented here will be directly for the 
infinite system. Details of the method can be found in 
ll|, [rj, |l3[ . In addition, we will als o ap ply Takahashi's 
modified spin-wave theory (MSWT) [lj| to this problem. 

In Fig. Q] the phase diagrams for the 5 = 1/2 isotropic 
AFM and FM Heisenberg models as defined in Eq. (Q]) 
are shown. The phase boundaries and order parame- 
ters are obtained using the TMRG algorithm. For the 
AFM we have a dimerized phase for any value of the 
elastic constant K/\J\ at low enough temperatures be- 
cause the gain in magnetic energy will always win. The 
phase transition is second order and the evolution of the 
order parameter is exemplified for K/\J\ = 2 in inset 
(b) of Fig. Q] For the FM, on the other hand, a dimer- 
ized phase exists only at finite temperatures and only 
if K/\ J\ < K C /\J\ ~ 0.118. Here we find a tricritical 
point (TCP) at (T TCP /| J|,X"tcp/| J1) - (0.696,0.116). 
For K < i^TCP the transition is first (second) order 
if T < T T cp (T > Ttcp), respectively. Inset (a) of 
Fig. Q] shows that the order parameter for K/\J\ = 0.1 
evolves indeed continuously at the upper phase boundary 
although it increases very steeply to one. Note that in 
the small window Ktcv < K < K c both the high and 



the low temperature transition will be first order. 

Next, we discuss the application of Takahashi's MSWT 
to this problem. Usual spin-wave theory is modified 
by introducing a Lagrange multiplier which enforces a 
nonmagnetic state at finite temperature. This guaran- 
tees that the Mermin-Wagner theorem is respected. For 
the isotropic FM chain, results obtained by MSWT have 
been shown to be in excellent agreement at low temper- 
atures with exact results obtained by the Bethe ansatz 
141 Il5| . For the dimerized chain the unit cell is doubled 
so that a Holstein-Primakoff transformation with differ- 
ent bosonic operators on the two sublattices is required. 
The diagonalized Hamiltonian in linear spin-wave theory 



is then given by H mag = Ne + Y,k { w k@kPk 
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with eo == JS 2 and the two magnon branches w\ 
2| J\S ( 1 ± v cos 2 k + S 2 sin 2 k ) . The constraint of zero 

magnetization NS = ~^2 k {nsi^k) + n s( w /f )} i s imple- 
mented by a Lagrange multiplier \x which acts as a chem- 
ical potential with ^s(^ fe ) = {exp[(cj* — fi)/T] — 1} _1 
being the Bose factors. For t/(l — 6 2 ) <C 1, where 
t = T/(\J\S) is the reduced temperature, we find an- 
alytically AS 2 ii/T = -t/(l - 5 2 ) + 0([t/(l - <5 2 )] 3 / 2 ). 
In the same limit the free energy per site is given by 
(/ - e )/T = a[t/(l - <5 2 )]!/2 + ([t/(l - S 2 )}), with 
a = -C(3/2)/(2v4r). For the FM chain we have there- 
fore a gain in magnetic energy due to a dimerization 
~ -T 3 / 2 <5 2 . 

To calculate spin correlation functions it is essential 
to take also quartic bosonic terms into account. For the 
bond correlations B s ^ = (S^'S^jii) this leads to 



b S ( W ) = (jj J2 {" s K ) - n B^i)} /k 



(2) 



with ff = (cos 2 k ± 5 sin 2 k)/\/cos 2 k + S 2 sin 2 k. The 
plus (minus) sign in (S , 2j-S , 2j±i) and in f k applies for 
the strong (weak) bond, respectively. We define 



Agg — (S2-J-S' 
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with A ss acting as an order parameter for the dimerized 
chain. In Fig. [U the MSWT and TMRG results for A ss 
are compared for the case of S = 1. The agreement is 
good for temperatures up to T/\ J\ ~ 1, in particular for 
small S. We also note that the MSWT gives a value in the 
fully dimerized case (5 = 1) which is in good agreement 
with the exact result, however, it predicts corrections for 
S = 1 — e (e -C 1) to be of order e 2 , whereas the numerical 
results and perturbation theory show that the corrections 
are of order e. In the inset of Fig. [5] it is shown that the 
phase diagrams for model ([I]) with 5 = 1/2 and 5 = 1 
are almost identical, if the axes are scaled appropriately. 

In Fig.[3]Ja) the correlation functions on the strong and 
weak bond for 5 = 1/2 are shown separately as a function 
of temperature for different d. We want to emphasize 




FIG. 2: (color online) Order parameter A ss (3[ as a func- 
tion of S for the S — 1 dimerized FM. The lines (symbols) 
denote the MSWT (TMRG) results, respectively. The inset 
shows the phase diagrams obtained by TMRG for the Hamil- 
tonian JT} with 5=1/2 (solid line) and S = 1 (dashed line) 
with both axes scaled appropriately. The phase transition is 
first (second) order for T < Ttcp (T > Ttcp)- The TCP is 
marked by a dot (square) for S — 1/2 (S — 1). 
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FIG. 3: (color online) TMRG results for S = 1/2: (a) B a(w) 
for the strong (solid line) and weak bond (dashed line), (b) 
A^ s for the same values of 5 as in (a), (c) the entanglement 
entropy SJaJa for different S as a function of temperature, and 
(d) S^uu) as a function of S for two temperatures. 



The concurrence for p — an entanglement measure com- 
monly used at zero temperature — can be expressed as 
C2 3 ,2 3 ±i{~p) = 2max{0, |(5+5y ±1 )|-|l/4+(5| J .5| j±1 )|}. 
It is zero for FM correlations. More interesting is the be- 
havior of the entanglement entropy, S °^ = —Trp In p. It 
is again zero for the fully polarized ground state which is 
a pure state. At finite temperature we have for a — s(w) 



s: at =(B a - 



For T -> 0, B 



' ^ In (- - B a ) - (B a 

(«,) - 1/4 and S^ - 



iJHl + T 

(4) 
In 3, see Fig. O^c). 



S° nt .\ therefore jumps signaling the phase transition at 



s(w) 

T = 0. For T 



oo, on the other hand, B 



s(w) 







and S^K — * 2 In 2. Quite generally, the entanglement 
entropy for a segment with n sites will go to nSx for 
T — ► oo, where St is the thermal entropy per site [1 SI ] - 
At any fixed finite temperature the entanglement entropy 
^stw) decreases (increases) on the strong (weak) bond 
with increasing modulation 5, see Fig. [3Jd). The gain in 
magnetic energy at finite temperature due to a dimeriza- 
tion might therefore also be seen as a gain in entangle- 
ment entropy on the weak bonds. 

Let us finally discuss the relevance of a thermally 
driven dimerization for systems with orbital degrees of 
freedom. This mechanism is particularly important for 
transition metal oxides with perovskite structure where 
the valence electrons are situated in the t^g orbitals. Be- 
cause tig orbitals are not bond oriented the electron- 
phonon coupling is weak so that we might ignore lattice 
degrees of freedom to first approximation. With appro- 
priately rescaled parameters, the physics discussed below 
is almost independent of the spin value S. For definite- 
ness, we will consider in the following the case of an ef- 
fective spin S = 1 appropriate for systems with a 3d 2 
valence electron configuration, as for example, YVO3, 
and a twofold orbital degeneracy described by an orbital 
pseudospin t = 1/2. A ID Harniltonian reflecting the 
spin-orbital physics for such a system is given by [191 ] 



H 3T ^jJ2(S r S j+1 + l)[r r T i+l 7 -•:•/; j. (•-,) 
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where J > is the superexchange and jh is proportional 
to the Hund's coupling and promotes FM spin correla- 
tions. Using a MF decoupling, which is reasonable for 
FM spin correlations [20|, we write H$ T — H$ + H T , 
where H$ (H T ) is the Harniltonian for the spin (orbital) 
sector, respectively. If we allow for a dimerization in 
goes through a maximum at some finite temperature, and both sectors then H S i T ) is — up to a constant — given 



again that for 6 £ [0,1) the ground state is still the usual 
FM state and the correlations on the weak and strong 
bond are thus identical, B s = B w = 1/4. The difference 
between the correlations on the strong and on the weak 
bond, A^ s , shown in Fig.[3jb) is therefore zero at T = 0, 



goes to zero again for T 



where B, 



(w) 



0. 



Another way of looking at the response of the FM chain 
to a periodic modulation is to study the entanglement of 
a weak or a strong bond with the rest of the system. Here 
we will concentrate on the case S = 1/2. The entries of 
the two-qubit reduced density matrix p for a bond can be 
related to the correlation functions on that bond [16|, [l7 1 . 



S(r) 

by Eq. (p} with J — ► JJs( T ), & ~ * $s{t)i an d S represent- 
ing the spins S — 1 or the orbital pseudospins r = 1/2, 
respectively. The effective superexchange constants are 
given by J s = A+ r /2 + 1/4 - lH and J T = A+ s /2 + 1, 
with A* r defined analogously to A ss . Strong quan- 
tum fluctuations for pseudospin 1/2 and jh > will 
favor AFM coupled orbitals, J T > 0, and FM coupled 
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FIG. 4: (color online) Phase diagram and dimerization pa- 
rameters Ss, S T for the spin-orbital model ([5]) with 7h =0.1 
in MF decoupling. The dashed lines denote the phase bound- 
aries between the uniform (U) and dimerized (D) phases. 



spins, Js < 0. The dimcrizations are then given by 
S s = A- T /{2J S ) and S T = A SS /{2J T ). This means that 
the exchange constants and the dimerizations for each 
sector are determined by the nearest-neighbor correla- 
tions in the other sector and therefore have to be calcu- 
lated self-consistently. We can simplify this procedure 
by noting that A"t„, , show only a weak dependence 
on dimerization and temperature for low temperatures. 
We therefore fix Js(t) by using the values for Aj„ , 
obtained for an undimerized chain at zero temperature. 
This leads to J T = 2 and J s = 1/2 - In 2 - j H [H|. 
Now the dimcrizations Ss( T ) can be easily determined 
self-consistently. The results for ju = 0.1 — which is a 
realistic value for cubic vanadates — are shown in Fig. [4] 
For 0.10 < T/J < 0.49 the self-consistent MF decoupling 



leads to nonzero values for S 



S(r)- 



The evolution of the 



dimerization parameters in this temperature regime has 
a dome-shaped form with a maximum at T/J ~ 0.2. In 
agreement with Fig. [TJ the dimerization in the FM spin 
chain is much larger than the dimerization in the AFM 
orbital chain and at T/J = 0.2 we have 5s ~ 0.86 which 
is already close to perfect dimerization (Fig. [4|. This 
underlines that the thermally activated dimerization in 
the FM chain is the driving force behind the finite tem- 
perature dimerized phase for the spin-orbital chain. The 
phase transitions at finite temperature between a uni- 
form and a dimerized phase are a consequence of the 
MF decoupling. Such phase transitions will not occur 
for the strictly ID model ©. Nevertheless, numerical 
calculations for this model [6j show that a dimerization 
is the leading instability at temperatures which support 
the dimerized phase in the MF decoupling solution. 

Summarizing, we have shown that a dimerization can 
occur in FM spin chains but has to be activated by ther- 
mal fluctuations. The gain in magnetic energy at finite 
temperatures can be related to an increased entangle- 
ment entropy on the weak bonds. For a FM chain with 
spin-phonon coupling we have derived the phase dia- 



grams as a function of temperature T and the effective 
elastic constant K for spin values S — 1/2 and 5 = 1. 
Thermodynamic properties of the dimerized FM chain 
can be calculated analytically with good accuracy for 
temperatures T < \J\S by a MSWT. Remarkably, this 
approach works for all dimcrizations 6 G [0, 1] if quartic 
terms are taken into account appropriately. For a system 
of coupled FM spin-1 and AFM orbital pseudospin-1/2 
degrees of freedom we found, using a mean-field decou- 
pling, a finite temperature dimerized phase. This shows 
that a dimerization is a universal instability of FM chains 
at finite temperatures, and may be triggered by the cou- 
pling to purely electronic degrees of freedom. This lat- 
ter mechanism seems to be relevant for many transition 
metal oxides with (nearly) degenerate orbital states. 
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